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PREFACE 


This  interim  report  was  prepared  by  Hughes  Research  Laboratories, 
Malibu,  California  under  contract  F30602 -76  -C-0021.  It  describes  work 
performed  from  1 August  1975  to  31  October  1975.  The  principal  scientist 
is  Dr.  James  E.  Pearson  and  the  principal  investigators  are 
Dr.  Melvin  E.  Pedinoff  and  Mr.  Stan  A.  Kokorowski.  This  project  is  a part 
of  the  Adaptive  Optics  program  in  the  Opto -Electronics  Department, 
managed  by  Dr.  Viktor  Evtuhov,  at  Hughes  Research  Laboratories. 


SUMMARY 


The  objectives  of  this  program  are  threefold:  (1)  to  determine 
what  kind  of  modulations  interfere  with  the  proper  operation  of  a multi  - 
dither  COAT  system  and  to  characterize  the  amplitude,  frequency  spec- 
trum and  spectral  power  density  of  these  modulations;  (2)  to  determine 
the  characteristics  of  receiver  modulations  produced  by  backscatter  from 
real  targets  in  expected  operational  scenarios  and  to  assess  their  effect- 
iveness in  degrading  COAT  system  convergence  and  performance;  (3)  to 
determine  an  optimum  COAT  system  design  that  is  insensitive  to  speckle 
induced  receiver  modulati  "H  to  define  the  critical  design  parameters. 
This  report  summarizes  the  w^  k performed  from  1 August  1975  to 
31  October  1975. 

We  have  experimentally  studied  the  speckle  patterns  generated 
by  scaled  realistic  target  model ) consisting  of  metallized  or  painted 
spheres,  cylinders,  and  cones  using  a helium -neon  (6328  A)  target- 
speckle  measurement  apparatus  specifically  assembled  for  this  program. 
Target  signatures  were  recorded  with  two  receiver  apertures,  one  equal  to 
the  transmitter  aperture,  and  one  equal  to  one -third  of  the  transmitter 
aperture.  The  larger  aperture  size  significantly  reduces  the  observed 
speckle  contrast  ratio  and  reduces  the  amplitude  of  the  higher  speckle 
spatial  frequency  components. 

Speckle  pattern  data  and  fast-Fourier-transform  spatial  frequency 
spectra  obtained  by  rotation  of  the  target,  or  by  translation  of  the  receiver, 
have  proven  to  be  identical  in  the  study,  thus  simplifying  the  remaining 
experimental  tasks  in  this  program. 

Target  speckle  signatures  taken  with  most  of  the  metallized  target 
shapes  are  very  similar,  with  the  exception  of  a cone-shaped  target 
illuminated  nose-oi  , This  target  showed  a very  broad  speckle  pattern  due 
to  the  small  scattering  cross  section  of  the  nose.  Experiments  with  white 
painted  (Lambertian  scat  ”ing)  targets  showed  reduced  contrast  ratios 
due  to  decorrelation  of  the  two  orthogonally  polarized  speckle  patterns. 
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The  madmum  contrast  ratio  observed  in  all  of  the  experiments  was  0.  89 
and  the  average  of  all  of  the  maximum  contrast  ratios  was  0.  56.  The 
average  contrast  ratios  observed  were  considerably  lower  and  are  estimated 
to  be  approximately  60%  of  these  figures. 

Real  time  spectrum  analyzer  spectra  obtained  by  rapid  rotation  of 
the  targets  are  identical  to  spatial  frequency  spectra  obtained  by  performing 
a fast  Fourier  transform  (FFT)  on  the  speckle  pattern  data.  Increasing  the 
target  rotation  speed  caused  the  speckle  spectrum  to  broaden  toward  high 
frequencies  in  agreement  with  the  FFT  spatial  frequency  results.  The 
speckle  signatures  produced  by  rotation  of  the  metallized  targets  are  identical 
with  those  produced  on  an  earlier  program  (contract  F30602-75-C-0001)  by 
rotating  scotchlite  targets.  Since  those  targets  failed  to  produce  severe 
degradation  of  COAT  system  performance,  the  targets  measured  and 
reported  here  are  expected  to  have  littfe  effect  on  system  performance. 

A statistical  analytical  model  o'.'  speckle  effects  in  multidither 
COAT  systems  has  been  developed.  The  model  has  shown  that  the  ac  signal 
power  generated  by  speckle  in  the  receiver  can  not  exceed  0.  5 times  the  dc 
power.  A speckle  amplitude  coefficient  Cb  has  been  defined,  equal  to  the 
square  root  of  the  ac  speckle  signal  power,  and  the  extreme  of  convergence 
as  a function  of  speckle  coefficient  has  been  calculated  for  comparison  with 
computer  simulation  results. 

Computer  simulations  of  multidither  COAT  performance  have  been 
made  using  an  ad  hoc  "speckle"  spectrum  and  using  theoretical  target  speckle 
data  generated  by  General  Research  Corporation.  Three  simulations  were 
run  corresponding  to  rotation  of  a uniformly  illuminated  spherical  target  at 
rates  of  0.01,  2,  and  10  rad  sec.  The  corresponding  speckle  coefficients  in 
our  analysis  are  0.  04,  0.  6,  and  0.  5.  The  latter  two  cases  are  sufficiently 
close  to  the  largest  possible  speckle  modulation  coefficient  of  0.  707  to  con- 
stitute worst  cases.  The  response  of  an  18 -channel  COAT  system  was 
partially  degraded  at  0.01  rad/sec,  and  severely  degraded  for  the  other  two 
cases.  Investigation  of  the  COAT  computer  model  revealed  significant 
second  harmonic  generation  in  the  AGC  at  65%  modulation  depth.  This  AGC 
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model  was  determined  to  be  unrealistic.  An  improved  AGC  model  in  the 
simulation  . led  to  good  ag.-eement  between  computer -predicted  convergence 
levels  and  'is  experimentally  observed  with  the  DARPA/RADC  COAT 
system.  The  simulation  model  is  not  finalized,  since  at  the  moment  it  pre- 
dicts longer  convergence  times  with  these  changes  than  without  them.  By 
comparing  the  statistical  analysis  to  the  simulation  results,  we  have  found 
that  the  statistical  analysis  serves  as  a good  predictor  of  ‘he  lower  level  of 
system  convergence. 

During  the  remainder  of  this  program,  experimental  measurements 
will  be  performed  with  the  DARPA/RADC  COAT  system  on  the  speckle- 
characterized  targets  under  conditions  of  thermal  blooming  and  turbulence. 

An  acousto -optic  speckle  simulation  experiment  will  provide  data  on  the 
minimum  modulation  index  required  to  disable  a single  COAT  servo  channel. 
Work  on  the  computer  simulation  model,  in  conjunction  with  the  General 
Research  Corporation,  will  continue  throughout  the  remainder  of  the  program. 
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I. 


INTRODUCTION 


A.  Program  Objectives 

The  overall  purpose  of  this  program  is  to  produce  a consistent, 
unified  theory  on  multidither  COAT/target  interactive  effects  that  correlate 
experimental  data  with  analytical  models.  Specifically  there  are  three  pri- 
mary objectives.  The  first  objective  is  to  determine  what  kind  of  modulations 
interfere  witn  the  proper  operation  of  a multidither  COAT  system.  The  signi- 
ficant parameters  of  the  modulation  are  amplitude,  frequency  spectrum,  and 
power  spectral  density.  The  second  objective  is  to  determine  the  character- 
ises of  receiver  modulations  pioduced  by  backscatter  from  real  targets  in 
expected  operational  scenarios;  that  is,  to  determine  which  ocenarios  can 
produce  the  modulations  which  will  affect  the  COAT  system.  The  third 
objective  is  to  determine  an  optimum  COAT  system  design  that  is  insensitive 
to  speckle -induced  receiver  modulations  and  to  define  the  critical  design 
parameters. 

B.  Research  Program  Plan 

The  18-element  multidither  COAT  system  developed  and  tested  on 
DAR PA/ R ADC  Contract  F30602 -73 -C-0248  was  expanded  to  a 21 -element 
system  to  include  autotracking  and  focus  controls  under  Amendment  No.  1 of 
DARPA/RADC  Contract  F30602-75-C-0001.  This  same  system  is  used  for 
laboratory  experiments  on  this  program,  which  is  in-depth  investigation 
of  speckle  phenomena. 

The  research  program  for  this  contract,  illustrated  in  Figure  1, 
runs  from  1 August  1975  to  1 February  1976.  Because  of  ongoing  measure- 
ments on  the  COAT  Measurements  and  Analysis  program,  a separate  target- 
speckle  measurement  apparatus  was  planned  including  a helium  neon  6328  A 
target  illumination  system. 
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The  program  is  divided  into  three  tasks:  Target  Backscatter 
Measurements,  COAT  System  Measurements,  and  Analytical  Studies.  Task  1 
represents  the  acquisition  and  characterization  of  targets  and  target  speckle 
data  not  previously  available  to  the  program.  Task  2 will  determine  COAT 
system  performance  in  the  presence  of  turbulence,  blooming,  and  speckle 
noise  modulation  from  these  targets.  Task  3 is  a computer  simulation  study 
designed  to  corroborate  or  predict  the  experimental  finding, 


TASK  1 TARGET  BACKSCATTER  MEASUREMENTS 

TASK  2 COAT  SYSTEM  MEASUREMENTS 

1.  INDUCED  MODULATIONS  IN  RECEIVER:  NO 

COAT  CORRECTION,  NO  DISTORTIONS 

2.  ACOUSTO-OPTIC  MODULATOR  EXPERIMENTS 


3.  COAT  CORRECTION,  SCALED  TARGET  SCENARIOS, 
TURBULENCE  AND  BLOOMING 

TASK  3 ANALYTICAL  STUDIES 


Figure  1.  COAT /Target-Signature  interaction  schedule 


II. 


SPECKLE  AND  TARGET  SIGNATURES 


The  physical  problem  known  as  "speckle  noise"  can  be  stated  briefly. 
When  a coherent  laser  beam  illuminates  a target,  some  of  the  energy  is 
scattered  back  in  the  direction  of  the  transmitted  beam.  The  backscattered 
radiation  produces  a random  intensity  pattern  (a  "speckle  pattern")  at  a 
receiver,  which  it  located  close  to  the  transmitter  for  cases  o'  interest  to 
us.  When  the  appropriate  conditions  exist,  any  movement  of  this  speckle 
pattern  relative  to  the  receiver  will  produce  an  amplitude  modulation  of  the 
received  signal.  Since  a multidither  COAT  system  receiver  senses  amplitude 
modulations  on  the  backscattered  beam  that  are  produced  by  dithering  the 
transmitted  beam  phase,  spurious  amplitude  modulations  in  the  receiver  may 
swamp  the  desired  modulations  and  thus  interfere  with  the  system  operation. 

Doppler  shifts  produced  by  rapidly  rotating  targets  may  also  produce 
false  signals,  but  these  spurious  signals  are  expected  to  be  less  important 
than  amplitude  modulation  effects  for  the  target  scenarios  of  interest.  r his 
is  true  because  the  most  likely  glint  point  on  a rapidly  rotating  target  is  one 
normal  to  the  beam  incidence  direction;  this  point  vill  have  nearly  zero 
transverse  doppler  shift.  Longitudinal  doppler  shifts  are  very  narrow  band 
and  thus  will  have  minimal  effect  on  a multiple -channel  COAT  system. 

The  speckle  characteristics  of  greatest  interest  to  this  program  are 
the  speckle  contrast  ratio  and  the  spatial  frequency  spectrum  of  the  speckle 
because  these  parameters  can  be  directly  translated  into  modulation  depth 
and  temporal  frequency  spectra  for  moving  targets.  These  parameters  are 
discussed  in  this  section  in  terms  of  both  analytical  models  and  experimental 

measurements. 
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A. 


Analytical  Speckle  Models 
1 . Models  for  Transmissive  Diffusers 

In  the  analytical  speckle  theory  derived  by  George,  et  al.  , a 
contrast  ratio,  CR,  is  defined  oy  the  ratio  of  the  standard  deviation  of  the 
speckle  intercity  to  the  average  intensity: 


where  the  expected  value  of  the  intensity  is  | ux  | and  cr(x)  is  its  standard 
deviation,  (In  the  analysis  of  the  experimental  results  presented  later  in  this 
report,  a heuristic  contrast  ratio  CR  , was  used 


C 


R 


1 


u 


xi 


u 


JLL 


u ■ + u . 
XI  Xj 


(2) 


where  u . and  u^j  are  the  instantaneous  adjacent  maximum  and  minimum 
values  of  a speckle  intensity  pattern.  This  datum  represented  the  maximum 

contrast  ratio  observed  in  each  experiment.) 

In  this  theory,  a transmissive  diffuser  with  total  surface  proportional 
to  A generates  a specular  transmission  scattering  term  A-BN  and  a roughness 
term  BN.  Thus  the  electric  field  received  at  X in  the  image  plane  is  given  by 

-in  h 

v(X)  = A-BN  + B ^ e m » (3) 

m-1 


where  h represents  the  random  height  within  the  element  and  h represents 
the  normalized  radian  frequency  variable  (phase  variable).  In  the  absence 


of  roughness  all  hm  = 0 and  v(X)  = A.  In  the  absence  of  specular  behavior, 
the  maximum  value  of  BN  is  A.  In  general,  the  specular  fraction  is  given  by 


S 


A - BN 


A 


(*> 


and  the  diffuse  fraction  is  given  by 


BN 

A 


(5) 


For  the  special  case  of  a very  rough  surface,  George's  expression  for  the 
contrast  x“atio  reduces  to 


(6) 


Thus  the  contrast  ratio  for  a very  rough  surface  depends  only  on  the  diffuse 
fraction,  b,  the  specular  fraction  S,  and  the  total  number  of  scatterers  N, 
and  not  on  the  detailed  density  function  for  height  fih^) 

For  the  nonspecular  case,  S = 0 (b  = 1)  and  the  contrast  is  1.  A 
more  general  calculation  of  the  contrast  ratio  as  a function  of  a frequency 
roughness  parameter  TK/tt  is  shown  in  Fig.  2.  The  variable  h is  the  root 
mean  square  height  of  the  diffuser  arid  the  parameter  >1  is  a normalized 
radian  frequency  variable.  In  effect,  H is  the  average  wave  number  differen- 
tial encountered  by  the  light  in  passing  through  the  phase  screen  (diffuser). 
Thus  the  product  Th/rr  for  the  random  phase  screen  is  identical  to  a surface 
roughness  parameter,  h/X.  for  reflection-generated  speckle. 
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ROUGHNESS  PARAMETER  OR  -~ 


Figure  2 


. Contrast  ratio  Cj^  versus  surface  roughness  parame.er. 
The  fractional  diffuse  reflectivity  component  is  b where 
s+b  = 1 (eq.  21)).  The  number  of  illuminated  scatterers 
is  m . 
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Several  important  conclusions  can  be  derived  from  Fig,  2,  First,  for 
a large  number  of  scatterers  (m  greater  tha..  - r equal  to  100),  the  specular 
component  need  only  be  increased  to  20%  (b  = 0.  8)  to  reduce  the  contrast 
ratio  to  0.  5.  Second,  for  fixed  values  of  b,  increasing  m almost  always 
decreases  the  contrast  ratio.  The  limited  exception  is  the  case  b = 1, 
h/X  > 0.  5 where  a slight  increase  in  contrast  ratio  occurs.  Finally,  for  a 
given  surface,  as  the  wavelength  increases,  the  contrast  ratio  decreac^s 
sharply  when  X exceeds  h.  These  results,  combined  with  the  fact  that 
surfaces  tend  to  become  more  specular  (smaller  values  of  b in  Fig.  2)  at 
longer  wavelengths,  indicate  that  speckle  modulation  effects  will  be  leas 
severe  at  10.  6 pm  or  3.  8 pm  than  at  0.  5 pm. 

2.  Probability  Density  Model  of  Speckle 

From  J.  M.  Burch,  ^ the  probability  density  function  P^fu^)  for 
a speckle  intensity  u^  containing  a spev  ular  component  is  given  by 


where  u^  is  the  mean  of  intensity,  is  the  contrast  ratio  defined  as  the 

ratio  of  the  standard  deviation  of  intensity  <r  to  the  mean  u , i.  e.  = °7UX. 

The  function  I is  the  modified  Bessel's  function  of  order  zero.  One  finds 
o 

that  for  C0  = 1,  the  density  function  reduces  to 
rv 


P (u  ) 

X X 


1 

— exp 
u 

X 


(8) 
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This  result  signifies  that  at  zero  intensity,  ux/ux  = 0,  and  for  high  contrast 
ratios  the  probability  density  of  zero  intensity  is  not  zero.  In  other  terms, 
in  the  absence  of  a specular  component  the  special  case  = 1 occurs,  and 
there  is  a high  probability  of  finding  an  intensity  zero  (null)  in  the  speckle 
pattern. 

For  unpolarized  or  depolarized  Lambertian  scattering,  the  speckle 

intensity  u is  the  sum  of  the  intensities  of  2 orthogonal  polarizations  ux  + uy 

which  are  uncorrelated.  The  density  PQ  for  the  total  intensity  is  calculated 

from  a convolution  of  P and  P . 

x y 

ru 

Po(u)  = J Px(u  - uy)  Py(uy)duy  . (9) 

o 

The  quantities  Pq(u)  and  Px(ux)  i*1  eqs.  (?)  and  (9)  are  plotted  in 
Figs.  3(a)  and  3(b).  The  contrast  ratios  used  to  produce  the  curves  in  Fig.  3 
were  taken  from  experimental  measurements  on  a series  of  diffusers  of  vary- 
ing thickness  by  George,  et  al.  1 A very  significant  conclusion  that  can  be 
drawn  from  Fig.  3(a)  for  a surface  roughness  corresponding  to  or  greater 
than  100  pm  diffuser  thickness,  is  a large  probability  density  at  zero 
intensity.  That  is,  the  polarized  specxle  patterns  will  have  deep  nulls  and 
contrast  ratios  of  1.  For  smaller  surface  roughness,  the  contrast  ratio 
will  be  less  than  1.  In  Fig.  3(b),  which  is  for  unpolarized  scatterings  (as 
from  a Lambertia.  surface),  contrast  ratios  of  1 are  found  only  for  a narrow 
range  of  the  roughness  corresponding  to  a diffuser  thickness  of  100  |*m. 

One  can  onclude  from  these  results  that  contrast  ratios  of  1 will  be 
observed  occasionally  with  polarized  receivers  when  the  surface  roughness 
is  large,  and  contrast  ratios  less  than  1 will  occur  with  unpolarized  receivers. 
This  conclusion  is  valid,  however,  only  for  targets  which  produce  significant 
depolarization  of  the  backscattered  energy.  In  visible  wavelength  experiments, 
Lambertian  scatterers  usually  exhibit  this  behavior.  Such  depolarization, 
however,  may  not  apply  at  3.  8 or  10.  6 (Jim  wavelengths.  We  must  hLa o be 
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This  result  signifies  that  at  zero  intensity,  ux/ux  = 0,  and  for  high  contrast 
ratios  the  probability  density  of  zero  intensity  is  not  zero.  In  other  terms, 
in  the  absence  of  a specular  component  the  special  case  Cn  = 1 occurs,  and 

X\ 

there  is  a high  orobability  of  finding  an  intensity  zero  (null)  in  the  speckle 
pattern. 

For  unpolarized  or  depolarized  Lambertian  scattering,  the  speckle 
intensity  u is  the  sum  of  the  intensities  of  2 orthogonal  polarizations  u^  + u 
which  are  uncorrelated.  The  density  PQ  for  the  total  intensity  is  calculated 
fiom  a convolution  of  P and  P . 

X y 

Po(u)  = J Px(u  - uy)  Py(uy)duy  . (9) 

o 

The  quantities  PQ(u)  and  Px(ux)  in  eqs.  (7)  and  '?)  are  plotted  in 
Figs.  3(a)  and  3(b).  The  contrast  ratios  used  to  produce  the  curves  in  Fig.  3 
were  taken  from  experimental  measurements  on  a series  of  diffusers  of  vary- 
ing thickness  by  George,  at  al.  * A very  significant  conclusion  that  can  be 
drawn  from  Fig.  3(a)  for  a surface  roughness  corresponding  to  or  greater 
than  100  pm  diffuser  thickness,  is  a large  probability  density  at  zero 
intensity.  That  is,  the  polarized  speckle  patterns  will  have  deep  nulls  and 
contrast  ratios  of  1.  For  smaller  surface  roughness,  the  contrast  ratio 
will  be  less  than  1.  In  Fig.  3(b),  which  is  for  unpolarized  scatterings  (as 
from  a Lambertian  surface),  contrast  ratios  of  1 are  found  only  for  a narrow 
range  of  the  roughness  corresponding  to  a diffuser  thickness  of  100  pm. 

One  can  conclude  from  these  results  that  contrast  ratios  of  1 will  be 
observed  occasionally  with  polarized  receivers  when  the  surface  roughness 
is  large,  and  contrast  ratios  less  than  1 will  occur  with  unpolarized  receivers. 
This  conclusion  is  valid,  however,  only  for  targets  which  produce  significant 
depolarization  of  the  backscattered  energy.  In  visible  wavelength  experiments, 
Lambertian  scatterers  usually  exhibit  this  behavior.  Such  depolarization, 
however,  may  not  apply  at  3.  8 or  10.  6 pm  wavelengths.  We  must  also  be 


(b)  UNPOLARIZEO 


Figure  3.  Theoretical  density  functions 
for  polarized  and  unpolarized 
speckle  intensity.  A,  F -25  pm 
diffuser;  B,  G - 50  pm  diffuser; 
C,  H - 100  pm  diffuser;  D,  J - 
200  pm;  E,  K - 500  pm. 
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careful  when  drawing  conclusions  from  data  based  on  transmissive  scatterers; 
reflective  targets  can  exhibit  quite  different  properties.  In  particular, 
depolarization  and  "diffuse"  scattering  are  not  necessarily  coupled  for 
reflective  targets. 

3.  Spatial  Frequency  Spectra  of  Speckle 

The  spatial  frequency  spectrum  of  speckle  is  proportions1  to 
the  autocorrelation  of  the  squared  amplitude  of  the  illumination  function. 

For  this  scenario,  i.  e.  , target  illumination  with  a Gaussian  beam,  the  auto- 
correlation is  also  Gaussian  but  of  twice  the  width.  It  has  been  shown  by 
3 

N.  George  that  the  spatial  frequency  distribution  (or  power  density)  of 
speckle  for  the  one  dimensions,!  case  is  given  by 


S(f  ) oc  A2(0)  6(i  ) \r  + (\r  )2 
' x'  ' ' ' x'  o ' o 


a<X, ro,fx) 


a(X,  r , f ) 
' o x' 


(10) 


where  * denotes  the  autocorrelation  with  respect  to  the  spatial  frequency 
variable  fx>  X is  the  optical  wavelength  s.nd  rQ  is  the  target  range.  In  the 
lefthand  term,  A2(0)  represents  the  total  power  in  the  illuminating  beam  and 
the  function  a(\,  rQ,  f^)  represents  the  amplitude  distribution  in  the  illuminat- 
ing beam.  In  the  present  investigation,  the  illumination  intensity  is  Gaussian 
of  the  form 

_ i xi 

I(X)  = e 7 a2  = | MX)  j 2 . (11) 


By  "diffuse"  we  mean  a target  that  completely  randomizes  the  phase  of  an 
optical  wave  incident  on  it,  not  necessarily  a wide-angle  or  Lambertian 
scatterer. 
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The  autocorrelation  function  oi  this  Gaussian  is  also  a Gaussian  of 


twice  the  width  given  by 


I(X)  * I(X) 


r>'  "s/tT  e 


(12) 


Substituting  the  transformation  variable  for  A one  obtains  for  the  cpectrum 

,2 

l\r  t f 
I o x 

2 


S(f  ) = 
' x 


6(f  )\r 
' x o 


\rQ  ( T^fn  e 


(13) 


where  cr  “ 0.42  times  the  full  width  at  half  maximum  of  the  Gaussian  beam. 
Using  the  experimental  values  <r  = 2.  8 x 10  ^ m,  X.  * 0.633  x 10  pm,  and 
rQ  = 1.8  m,  we  obtain  the  amplitude  and  power  spectrum  results  shown  in 
Fig.  4.  When  finite  aperture  averaging  from  a 0.  15  mm  aperture  is  used  to 
weight  the  speckle  spectrum,  the  lower  curves  result.  We  will  see  later  that 
these  spatial  frequency  spectra  are  fairly  accurate  representations  of 
measured  phenomena. 

B.  Experimental  Speckle  Measurements 

Experimental  measurements  at  Caltech  are  shown  in  Support  of  the 
theory  by  George  et  al.  1 in  Figs.  5(a)  and  5(b).  The  general  features  of  the 
theoretical  curves  in  Fig.  3 and  the  experimental  curves  in  Fig.  5 are  in 
excellent  agreement.  The  polarized  receiver  case  with  the  500  pm  diffuser 
shows  a contrast  ratio  of  1 as  expected.  A simple  summary  of  the  thick 
diffuse  experimental  speckle  measurements  is  presented  in  Fig.  6.  For  our 
purpose,  the  diffuser  thickness  can  be  likened  to  surface  roughness.  lor  the 
unpolarized  receiver  case,  the  contrast  ratio  is  always  less  than  1 and  drops 
to  0.  5 for  large  thickness  (roughness).  The  maximum  contrast  ratio  of  0.  9 
occurs  at  a diffuser  thickness  of  100  pm  and  falls  off  on  either  side  for  the 
unpolarized  case.  For  reflective  targets,  this  is  equivalent  to  the  statement 
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SPECTRAL  DENSITY,  ARBITRARY  UNITS 


Figure  4.  Theoretical  spatial  frequency  spectra  of  speckle  with  Gaussian 
illumination,  (a)  Amplitude  spectrum  with  point  detector, 

(b)  Amplitude  spectrum  with  0.15  mm  detector,  (c)  Power 
spectrum  with  point  detector,  and  (d)  Power  spectrum  with 
0.15  mm  detector. 
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III. 


TFCHNICAL  ACCOMPLISHMENTS 


A.  Experimental  Target  Reflection  Speckle  Measurements 

1.  Preliminary  Experimental  Results 

During  the  latter  portion  of  the  previous  program  entitled, 
"COAT  Measurements  and  Analysis"  (RADC  F30602 -75-C-0001),  some  pre- 
liminary experimental  speckle  measurements  were  made,  motivated  by  the 
work  of  Ogroanik  and  Gurski.  4 This  work  is  reported  in  Quarterly  Technical 
Report  No.  3 for  the  period  of  1 January  1975  to  31  March  1975.  The  relevant 
results  obtained  in  that  study  were  that  a rotating  scotchlite  target  produced 
a significant  increase  in  COAT  receiver  speckle  noise  but  no  significant 
degradation  in  converged  power,  and  that  an  artificially  simulated  swept 
narrow  band  speckle  noise  spectrum  could  partially  degrade  COAT  system 
convergence. 

2.  Quasi  Static  Speckle  Measurements 

An  apparatus  shown  in  Fig.  7 was  assembled  to  measure  the 
speckle  characteristics  of  various  targets  developed  for  this  program.  A 
helium  neon  laser  beam  was  spatially  filtered,  truncated,  and  recollimated 
to  0.  15  mm  diameter  (i.  e.  , intensity)  to  simulate  the  converged  output  of  the 
COAT  system.  This  heam  then  illuminated  a variety  of  targets  at  a range  of 
1.8  m where  the  beam  radius  was  10  mm  to  the  first  null.  The  targets  were 
mounted  on  a motor-driven  tunntable  and  the  photomultiplier  detector  on  a 
motor  driven  translating  stage.  The  photomultiplier  (PMT) 
output  was  processed  in  a lock-in  amplifier  and  recorded  via 
a strip  chart  recorder.  A 100  Hz  synchronous  chopper  placed  in  the  path 
of  the  transmitted  laser  beam  provided  the  optical  modulation  and  the 
reference  signal  for  the  synchronous  detector.  In  the  series  of  experiments 
performed,  either  the  target  was  rotated  and  the  PMT  was  fixed,  or  the  PMT 
was  translated  and  the  target  was  held  fixed.  Another  parameter  that  was 
varied  during  the  experiment  was  the  size  of  the  receiving  aperture.  Pinholes 
of  either  0.  05  mm  or  0.  15  mm  were  placed  in  front  of  the  PMT  for  this 
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purpose.  The  larger  diameter  pinhole  corresponded  to  the  identical 
transmitter  and  receiver  aperture  diameter.  Data  taken  with  the  smaller 
aperture  showed  the  presence  of  higher  spatial  frequency  speckle  modulations 
which  are  integrated  out  by  the  larger  receiver  aperture.  The  PMT  receiver 
was  also  provided  with  a removable  polarizer  to  facilitate  measurements  of 

speckle  contrast  as  a function  of  polarization. 

Measurements  were  made  on  target  spheres,  cylinders,  and  cones 
fabricated  from  ordinary  glass.  Some  of  these  targets  were  sandblasted,  and 
others  were  left  smooth.  All  were  electroplated  with  nickel.  Some  of  these 
targets  were  then  coated  with  3M  NEXTEL  white  paint  which  shows  excellent 
Lambertian  scattering  properties. 


a.  Spatial  (Speckle)  Modulation  Generated  by  Target 

Rotation  - The  spatial  modulations  or  optical  intensity  patterns  produced  by 
slowly  rotating  metallic  spheres  and  cylinders  are  shown  in  Figs.  8,  9,  and 
10.  The  data  in  the  upper  curve  of  each  figure  was  taken  with  a 0.  15  mm  pin- 
hole aperture  receiver  and  in  the  lower  curve  with  a 0.  05  mm  aperture. 
Aperture -averaging  with  the  larger  receiver  is  clearly  evidenced  by  the  loss 
of  high  frequency  information.  The  data  taken  with  the  smaller  receiver 
(0.  05  mm)  shows  higher  spatial  frequency  content  and  a higher  contrast 

ratio  than  that  taken  with  the  larger  receiver  (0.  15  mm). 

The  maximum  contrast  ratios  for  the  smooth  sphere,  the  sandblasted 
sphere  and  the  smooth  cylinder  obtained  with  the  larger  aperture  were  0.  79, 
0.75,  and  0.64,  respectively,  and  0.  77,  0.  89  and  0.  80  when  measured  with 
the  smaller  aperture.  The  average  contrast  ratios  were  0.  73  and  0.  82  for 
the  large  and  small  aperture.  These  contrast  ratios  were  calculated  by  taking 
the  maximum  irradiance  and  an  adjacent  minimum  and  finding  the  ratio  of  the 
difference  to  the  sum.  The  average  contrast  ratios  observed  were  consider- 
ably lower.  Maximum  contrast  ratios  of  1. 0 corresponding  to  zero  nulls  were 
not  observed,  possibly  because  of  the  finite  receiver  size.  In  our  experiments, 
decreasing  the  aperture  size  increased  the  contrast  ratio;  therefore,  as  the 
aperture  becomes  small  compared  with  the  speckle  size,  the  maximum 
observed  contrast  may  approach  1. 
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Figure  9.  Speckle  patterns  generated  by  rotation  of  rough 
metallized  sphere. 

(a)  0.15  mm  receiver 

(b)  0.05  mm  receiver 


Figure  10.  Smooth  metallized  sphere  speckle  patterns 
generated  by  rotation 

(a)  0.15  mm  receiver 

( b)  0.05  mm  receiver 

After  a careful  study  of  these  data,  we  conclude  that  *he  shape  of  the 
target  and  the  nature  of  the  surface  roughness  does  not  significantly  affect 
the  spatial  modulation  spectrum  or  the  average  contrast  ratio,  as  long  as  the 
entire  illuminating  beam  is  smaller  than  the  target. 

b.  Equivalence  of  Rotational  and  Translational  Speckle 

Modulations  and  Spatial  Frequency  Spectra  — It  was  necessary  to  measure 
the  relationship  between  rotationally  and  translationally  generated  speckle 
modulations  to  simplify  the  number  of  COAT  convergence  experiments  re- 
quired in  the  next  phase  of  the  program.  The  translation  experiment  is 


particularly  troublesome,  since  it  cannot  be  performed  continuously  at 
constant  velocity  and  since  it  requir  3 twice  as  many  measurements  as  the 
rotational  experiment. 

We  have  measured  speckle  signatures  independently  by  detector 
translation  and  by  target  rotation.  These  two  measurements  are  not  equiva- 
lent because  rotation  of  the  target  through  a small  angle  produces  a small 
change  in  the  scattering  surface  near  the  periphery  of  the  beam.  As  long  as 
the  rotation  angles  are  small,  however,  this  is  a small  effect.  Several  data 
records,  in  fact,  exhibited  speckle  patterns  that  were  almost  mirror  images 
of  each  other  when  the  trans lation- rotation  functions  are  appropriately 
sequenced.  Figure  11  shows  an  example  of  this  type  of  data.  This  result 
establishes  that  for  small  angles  of  rotation  or  for  small  translations  two 
motions  were  equivalent. 

3.  Spatial  Frequency  Analysis 

The  spatial  frequency  spectrum  of  the  speckle  data  for  the 
smooth  sphere  shown  in  Fig.  8 has  been  calculated  via  a fast  Fourier 
transform  (FFT)  routine.  Figure  12(a)  represents  the  speckle  spectrum  of  a 
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Figure  11.  Symmetry  of  rotational  and  translational  speckle 
modulation  patterns. 
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Figure  12(b).  Spatial  frequency  spectrum  of 
smooth  metallized  sphere  — 
0.05  mm  receiver. 
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rotating  target  seen  by  a 0. 15  mm  receiver.  When  the  receiver  is 
translated,  the  general  shape  of  the  spectrum  is  very  similar  to  that  shown 
in  Fig.  12(b).  The  notch  in  the  spectrum  centered  about  0.67  x 10*M“*  is 
attributed  to  the  spatial  integration  of  the  signature  by  the  finite  aperture. 
Figure  12(b)  represents  the  same  spectrum  received  by  a 0.05  mm  receiver. 
Notice  the  improvement  of  the  higher  frequency  components  relative  to  the 
low'i . A naive  approach  to  the  effect  of  finite  aperture  integration  eliminates 
all  spatial  frequencies  higher  than  that  corresponding  to  the  receiver  aper- 
ture. The  calculated  data  for  a one  dimensional  slit  of  0. 15  mm  width  shown 
in  Fig.  12(a)  clearly  shows  the  presence  of  higher  spatial  frequencies  at  a 
level  13  dB  below  the  maximum.  Figure  12(b)  shows  the  effect  of  aperture 
integration  with  a receiver  one  third  the  size  of  the  transmitter  (0.05  mm). 
The  relative  effect  of  aperture  integration  is  that  the  larger  aperture  reduces 
the  high  frequency  spectrum  relative  to  the  low  by  9 dB. 

The  absolute  effect  of  aperture  integration  by  the  larger  aperture 
relative  to  a point  receiver  is  to  reduce  the  target  signature  spectrum  by  1 1 
to  13  dB  for  all  frequencies  beyond  the  aperture  cutoff  frequency.  Also,  all 
frequencies  greater  than  one-half  the  large  aperture  cutoff  frequency  are 
attenuated  by  more  than  3 dB . These  results  substantiate  our  experimental 
findings  in  Figs.  8,  9,  and  10  of  lowered  speckle  signature  contrast  when 
integration  by  the  larger  aperture  is  present. 

Signatures  of  sandblasted  metallized  spheres  with  an  integrating 
aperture  of  0. 15  mm  are  very  similar  to  the  smooth  sphere  case.  The 
spectrum  for  this  case  is  shown  in  Fig.  13  and  is  very  similar  to  those 
analyzed  for  the  smooth  sphere.  This  result  shows  that  the  nature  of  the 
speckle  spectrum  is  primarily  determined  by  the  target  illumination  function 
and  the  surface  shape,  and  not  by  the  state  of  the  surface  for  this  class  of 
specular  metallized  targets. 

Figures  14(a)  and  (b)  show  the  speckle  target  signatures  obtained  by 
rotating  a smooth  metallized  cone  about  an  axis  parallel  to  the  base  of  the 
cone.  This  signnature  is  the  anomally  of  the  data  taken  in  this  program. 

The  low  apparent  periodicity  is  partially  explained  by  the  geometry  of  the 
conical  surface.  The  only  significant  backscatter  appears  to  be  generated 
at  the  nose  of  the  cone.  This  implies  that  the  effective  dimension  of  the 
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Figure  14.  Speckle  pattern  from  smooth  metallized  cone  with  axial 
illumination. 

speckle  generator  is  smaller  than  the  illumination  function,  and  that  the 
speckle  spatial  period  should  be  large.  However,  the  speckle  pattern  is  not 
sinusoidal  but  highly  triangular  and  therefore  rich  in  spatial  frequency  har- 
monics. The  FFT  spectrum  is  shown  in  Fig.  15.  This  spectrum  drops 
from  the  dc  level  by  approximately  16  dB  at  high  frequencies.  The 
aperiodic  nature  of  the  spatial  function  explains  presence  of  high  frequen- 
cies in  the  spectrum.  The  linear  slope  of  the  spectrum  plotted  on  log  paper 
corresponds  to  an  exponential  or  Gaussian  frequency  discribut ion.  Thus, 
the  observed  spatial  frequency  spectrum  is  qualitatively  in  agreement  with 
the  theory  (see  Section  II-A,  eqs.  10  to  13). 

Target  velocity  and  rotation  rate  govern  the  real  time  frequency 
spectrum  generated  in  a receiver  by  speckle  modulation.  1 he  real  time 
spectrum  G(f  ) can  be  found  from  the  spatial  frequency  spectrum  G(F  ) by 
the  relationships 

ft  = <Fx>  • (Vx) 

o r 

ft  = <FX>  • <2<“V  • Z 
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where  V is  the  relative  velocity  of  the  receiver  and  the  target,  is  the 

relative  angular  velocity  of  the  target,  Z is  the  target  to  receiver  distance, 

f is  the  temporal  frequency  (sec"^)  and  F is  the  speckle  spatial  frequency 
1 1 X 

The  pertinent  parameters  for  our  scaled  laboratory  version  of  a 

typical  target  scenario  are  shown  in  Table  1.  As  an  example,  at  the 

4 -1 

experimental  spatial  frequency  midband  (F^  = 0.  5 x 10  m ).  the  temporal 
frequencies  ft  are  4 kHz  and  540  Hz  for  translation  and  rotation, 
respectively. 


Table  1.  Parameters  Determining  Speckle  Modulation  Spectrum 


Parameters 

Real  Scenario 

— 
Scaled  Laboratory 
Model 

Wavelength 

X - 3.  8 pm 

X = 0.  6328  pm 

Range 

Z = 4 x 103  m 

z = 1.8  m 

Transmitter  Aperture 

Drp  = 0.  7 m 

_ 

dt  = 0.  15x10  ' m 

Target  Translation 
Velocity 

V = 4 x 103 
8 / 

m/  sec 

v = V D /D„ 
s s t'  T 

= 8 x 10  * m/sec 

Target  Angular 
Velocity 

= 0.  06  rad/ sec 

X 

= (zDt/zDT).£!x 
= 0.  03  rad/sec 

4.  Target  Surface  Phenomena 

Figure  16  compares  the  speckle  modulation  signatures 
generated  by  rotating  a smooth  metallized  sphere  and  a white  painted  sphere. 
The  integrating  aperture  size  for  these  data  and  for  the  remainder  of  this 
section  is  0.  15  mm.  The  maximum  speckle  contrast,  ratios  for  the  smooth 
metallized  sphere  and  the  white  sphere  were  found  to  be  0.  79  and  0.  42, 
respectively.  The  Lambertian  scattering  characteristics  of  the  white  surface 
(known  from  measurements  on  other  programs)  was  expected  to  produce  less 
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specular  scattering  and  a higher  contrast  ratio.  However,  approximately  half 
of  the  scattered  energy  is  depolarized  and  the  speckle  pattern  for  the  orthog- 
onal polarization  is  uncorrelated  with  that  of  the  incident  polarization.  Fig- 
ure 16(c)  shows  the  effect  of  placing  a polarizer  in  front  of  the  detector 
oriented  parallel  to  the  transmitted  polarization.  The  maximum  contrast 
ratio  increased  to  0.  71.  This  result  qualitatively  substantiates  the  theory 
and  allows  us  to  predict  that  with  true  Lambertian  scatterers,  the  speckle 
contrast  ratio  will  be  substantially  less  than  0.  71,  provided  the,  receiver  is 
not  polarized  and  aperture  averaging  is  employed. 

One  of  the  more  interesting  target  shapes  studied  has  been  the  cone., 
The  results  of  a series  of  measurements  on  a cone  illuminated  from  the  side 
are  shown  in  Fig.  17.  In  the  case  of  the  smooth  metallized  cone  (Fig.  17(a)), 
the  low  maximum  contrast  ratio  of  0.  17  was  attributed  to  the  specular  nature 
of  the  cone  reflectance.  The  rougher  metallized  cone  surface  data  shown  in 
Fig.  17(b)  showed  a higher  contrast  ratio  of  0.43.  A contrast  ratio  of  0.58 
was  generated  by  the  white  cone  data  shown  in  Fig.  17(c). 

A similar  set  of  experimental  results  is  presented  in  Fig.  18  for  a 
cone-shaped  target  illuminated  head-on.  As  mentioned  before,  the  anoma- 
lously low  periodicity  in  Fig.  18(a)  is  explained  by  the  small  nose  cross 
section.  The  maximum  contrast  ratio  of  0.  79  is  consisted  with  a reflector 
which  produces  very  little  specular  scattering  in  the  direction  of  the  receiver 
because  of  its  geometrical  shape.  The  rough  cone  result  in  Fig.  18(b)  shows 
higher  spatial  frequency  content  and  lower  contrast  ratio  of  0.  56.  In 
Fig.  18(c),  the  white  cone  shows  similar  spatial  frequency  content  to  the 
rough  cone  of  Fig.  18(b),  except  that  the  maximum  contrast  ratio  of  0.35  is 
significantly  wer.  The  mechanism  of  Lambertian  depolarization  of  the 
target  return  at  normal  incidence  can  explain  this  contrast  reduction;  when 
these  targets  are  illuminated  with  the  laser  and  observed  through  a polarizer 
only  slight  depolarization  is  seen  for  the  metallized  surfaces,  but  total 
depolarization  is  observed  for  the  white  surfaces. 

A tabulation  of  all  of  the  measured  maximum  contrast  ratios  is  given 
in  Table  2.  The  largest  contrast  observed  with  equal  transmitter  and 
receiver  apertures  was  0.  79  and  the  average  of  the  maximum  was  0.  56.  The 
average  speckle  contrast  ratio  observed  in  any  experiment  was  of  the  order 
of  0.  5 or  less. 


DISTANCE.  nr»m 

Figure  18.  Speckle  patterns  of  cones  with  axial  illumination  and  nonaxial 
rotation 

(a)  Smooth  metallized  cone,  0.15  mm  aperture 

(b)  Rough  metallized  cone,  0 . 1 5 mm  ape rture 

(c)  White  painted  cone,  0.15  mm  aperture 


44 


i 


1 


5.  Dynamic  Speckle  Measurements 

The  real  time  amplitude  spectral  density  of  some  of  the  real- 
istic target  models  described  in  the  previous  section  were  measured  by 
rotating  the  targets  at  various  speeds.  The  receiver  signal  was  then  observed 
with  a spectrum  analyzer  (see  Pig.  7).  The  spectra  observed  generally 
agreed  with  the  preliminary  experimental  results  obtained  on  the  previous 

Table  2.  Maximum  Speckle  Contrast  Ratios  of  Variously 
Shaped  and  Surfaced  Targets 


Target 


Aperture, 

0.  15  mm 

Smooth  Sphere 

0.  79 

Rough  Sphere 

0.  75 

Smooth  Cylinder 

0.  64 

Average 

0.  73 

Maximum  Contrast  Ratio 


Aperture, 
0.  05  mm 


0.  77 
0.  89 
0.  80 

0.  82 


Smooth  Sphere 

White  Sphere 

White  Sphere  (Polarized) 


0.  79 
0.  42 
0.  71 


Smooth  Cone  (Side)  Illumination 
Rough  Cone  (Side)  Illumination 
White  Cone  (Side)  Illumination 


0.  17 
0.43 
0.  58 


Smooth  Cone  (End)  Illumination 
Rough  Cone  (End)  Illumination 
White  Cone  (End)  Illumination 

Average 


0.  79 
0.  56 
0.  35 

0.  56 
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contract^  with  a rotating  scotchlite  cylinder,  and  with  the  detailed  spatial 
frequency  spectra  extracted  from  the  realistic  target  speckle  signature  data. 

Figure  19  shows  the  spectrum  of  a smooth  metallized  spherical  target 
in  a stationary  condition  and  rotated  at  0.  67  revolutions  pel  second,  the  data 
having  been  taken  at  10  and  200  Hz  bandwidth.  Aside  from  a 0.  4 gain  factor 
for  the  200  Hz  bandwidth  case,  the  spectral  results  are  very  similar.  The 
noise  spectral  amplitude  due  to  speckle  is  approximately  a factor  of  2.8  larger 
than  the  photodetector  and  laser  noise  for  both  cases.  This  result  is  very 
similar  to  that  obtained  earlier  for  scotchlite  targets. 

A similar  set  of  test  results  is  shown  in  Fig.  20  for  the  smooth  metal- 
lized cvlindrical  target.  The  cylindrical  geometry  tends  to  produce  spread- 
ing cl  the  scattered  energy  in  only  one  coordinate  direction  as  compared  to 
th$  uniform  angular  scattering  from  the  spherical  target.  This  leads  to  a 
stronger  speckle  signal  as  evidenced  by  the  speckle  noise  being  a factor  of 
3.8  greater  than  the  system  noise. 

The  effect  of  varying  the  target  rotation  rate  is  shown  in  Figs.  21  and 
22  for  the  spherical  and  cylindrical  targets.  A definite  shift  of  the  amplitude 
spectrum  toward  high  frequencies  with  increasing  target  speed  is  evident  in 
both  cases.  The  larger  receiver  speckle  signal-to-no; se  ratio  produced  by 
the  cylinder  makes  the  effect  more  pronounced.  The  amplitude  spectral 
density  is  always  larger  at  the  low  frequency  end  but  increases  at  higher 
freq  mcies  with  increasing  target  rotation  speed.  In  other  words,  the  fre- 
quency of  the  3 dB  speckle  noise  level  shifts  proportionately  with  rotztion 
rate.  This  is  in  exact  agreement  with  the  discrete  spatial  frequency  spectral 
data  analysis  of  the  previous  section,  where  the  spatial  spectrum  is  always 
stronger  at  low  frequencies  and  falls  off  at  high  frequencies  due  to  the 
Gaussian  illumination  function  and  due  to  the  aperture  averaging  effect. 

Since  the  temporal  frequency  is  obtained  from  the  spatial  frequency  by 
multiplication  by  the  product  of  the  angular  rate  and  the  range,  this  spectrum 
is  either  expanded  or  compressed  as  the  velocity  of  the  target  is  changed, 
but  its  shape  is  basically  unchanged.  If  the  target  is  rotated  fast  enough,  the 
speckle  noise  spectrum  will  appear  flat  because  the  3 dB  cutoff  will  fall  outside 
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Figure  19.  Receiver  signal  amplitude  spectra  for  smooth  spherical 
metallized  target,  (a;  and  (c)  Stationary  target, 

(b)  and  (d)  Rotating  target.  The  spectrum  analyzer 
bandwidth  is  10  Hz  in  (a)  and  (b)  and  is  200  Hz  in 

(c)  and  (d).  Target  rotation  rate  0.67  rps. 


CTRAL  AMPLITUDE  SPECTRAL  AMPLITUDE 


FREQUENCY.  kHz  FREQUENCY.  kHz 

Receiver  signal  amplitude  spectra  for  smooth  cylindrical 
metallized  target.  <a)  and  {c)  Stationary  target,  (b)  and 
(d)  Rotating  target.  The  spectrum  analyzer  bandwidth  is 
10  Hz  in  (a)  and  (b)  and  is  200  Hz  in  (c)  and  (d).  Target 
rotation  rate  is  0.67  rps. 
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Figure  21.  Receiver  signal 
amplitude  spectra  for  smooth 
spherical  rotating  target. 

The  rotation  rates  are 
(a)  0.11  rps,  (b)  0.25  rps, 

(c)  0.38  rps  . The 
receiver  bandwidth  is 
200  Hz. 
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Figure  22.  Receiver  signal 
amplitude  spectra  for 
smooth  cylindrical  metal- 
lized rotating  target. 

(a)  0.21  rps , (b)  0.33  rps, 
(c)  0.5  rps.  The  receiver 
bandwidth  is  200  Hz. 


the  ai.pl.y  d :r.c.,.l  range.  A table  of  .pecHe  rnodula.ion  frequencie.  aa  a 

function  of  targe,  rotation  rate.  i.  given  in  Table  i.  A dtvtding  line 
the  -3  dB  signal  amplitude  regime. 

Table  3.  Modulation  Frequency  ac  a Function  of 
xaDie  « _ .n>)  nrttaf-inn  Rate 


1 

Rotation  Rate  -«*■ 
Spatial  Frequency 

0.25  RPS  0.5  RPS 

Modulation  Fret 

0.75  RPS  7 
juencies,  kHz 

1 RPS 

4 -1 

0x10  M 

0 

0 

0 

0 

4 -1 

0.25  x 10  M 

14 

28 

42 

' 

57 

4 ..-I 

0.50  x 10  M 

28 

57 

84 

113 

4 -1 

0.75  x 10  M 

42 

85 

127 

170 

j 

1.00  x 104  M"1 

57 

113 

170 

226 

1 

NOTE:  (Doublcline 

is  approximate  location  of 

-3  dB  amplitude  response.)  j 

A similar  se,  of  experiments  was  performed  with  a painted  sphertcal 
target,  known  to  be  a good  Lambertian  scatterer.  This  targe,  ytelded  ^trl 
nerg  in  both  return  polarisations  and  a very  low  speck.,  no.se  contr.bm 
when  recorded  with  an  unpolarised  receiver.  This  result  is  .dent.cal 
measured  in  the  quasi  static  case. 
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IV. 


SPECKLE  EFFECTS  ON  MULTIDITHER  SYSTEMS 


A.  Analytical  Model  of  Speckle  Modulation 

The  statistical  analysis  reported  earlier  has  been  generalized  to  be 
applicable  to  any  generalized  spsck’e  modulation  function,  Originally, 

the  whole  analysis  was  based  on  a speckle  modulation  function  with  a randomly 
chosen,  discrete  frequency  spectrum  centered  around  each  dither  frequency 
over  a very  narrow  bandwidth.  The  amplitude  of  this  modulation  was 
characterized  by  a parameter  a , and  COAT  system  convergence  levels 
were  predicted  as  a function  of  a . This  made  the  analytical  data  very 
easy  to  compare  to  computer  simulation  data.  With  a more  general  type  of 
modulation,  however,  the  parameter  a may  be  difficult  or  perhaps 

impossible  to  define.  Thus,  we  consider  the  very  general  type  of  speckle 
modulation  function  defined  by 


/co 

dw  a(o)  e‘iu,t  . (14) 

- co 

The  Fourier  transform  relationship  is 

a(w)  = — - J*  dtMg(t)oiu£  . (15) 

- CO 

To  be  physically  meaningful  M (t)  must  always  be  ru  ative. 

s 

Negative  values  imply  negative  intensity  a*  . receiver,  whicu  is  obviously 

nonphysical.  Conservation  of  energy  also  ’quires  that  the  average  value  of 

M (t)  be  unity.  We  illustrate  this  point  as  follows.  In  the  computer  Simula- 
s 

tions  of  the  COAT  servo  (see  next  section),  we  do  not  include  the  effect  of 
atmospheric  propagation.  The  receiver  signal  is  therefore  just  a constant 
multiplied  by  the  instantaneous  intensity  of  the  light  on  the  target  glint.  The 
instantaneous  energy  incident  on  the  receiver  is  thus  directly  proportional  to 


Preceding  page  Wank 


the  energy  incident  on  the  glint.  When  speckle  effects  are  introduced  into  the 
simulation  and  into  the  statistical  analysis,  we  simply  multiply  the  receiver 
signal  by  a modulation  function,  M (t).  Since  speckle  is  a v/ave  — interference 
effect,  as  opposed  to  an  absorption  or  scattering  effect,  it  must  be  unrelated 
to  the  total  energy  in  the  incident  or  reflected  fields.  In  other  words,  since 
speckle  is  an  energy  conserving  phenomena,  it  must  be  introduced  in  such  a 
way  that  the  total  energy  in  the  reflected  field  is  unaltered.  This  is  accom- 
plished in  the  computer  simulations  and  in  the  statistical  analysis  by  requiring 
that  the  average  energy  incident  on  the  receiver  remain  the  same  with  or 
without  the  mr  dulation  function  Mg(t).  The  intensity  at  the  receiver  is  no 
longer  required  to  be  instantaneously  proportional  to  the  intensity  on  the  tar- 
get glint  but  the  average  intensity  at  the  receiver  must  be  proportional  to 
the  average  intensity  on  the  glint,  and  the  constant  of  proportionality,  K, 
muft  still  be  the  same.  Mathematically  this  implies  that  the  mean  value  or 

dc  level  of  M (t)  must  be  unity, 
s 

Since  we  are  interested  in  the  power  ratio  of  the  COAT  system  dither 
signal  to  the  speckle  induced  signal,  we  first  calculate  the  average  power  flux 
at  the  receiver  due  to  speckle  modulations  alone,  (i.  e.  We  assume  the  usual 
"no-speckle"  intensity  is  unity).  While  the  total  dc  level  must  be  unity,  as 
discussed  above,  the  ac  power  at  the  dither  frequencies  is  yet  to  be  deter- 
mined. By  definition,  the  total  average  power  in  the  speckle  modulation, 


TS 


lim  1 
T— ► « 2T 


r 

/ 


Ms(t)  Ms(t)  dt 


(16) 


Using  the  definition  of  M (t)  in  eq.  (14)  gives: 

s 


> - lim 
TS  T-*>® 


a(u)  a (w)  dw 


(17) 
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A limiting  value  for  P occurs  when  M (t)  consists  of  a unit  dc  level 

Tb  8 

and  a unit  amplitude  single  frequency  modulation.  In  this  case, 


Ma(t)  = M(t) 


1 + sin  (wt  + w) 


(18) 


Carrying  out  the  integral  in  eq.  (16)  for  this  modulation  gives  Pg^  = 1.5. 

A more  realistic  modulation  function  would  include  several  frequencies 
with  random  phases.  By  increasing  the  number  of  frequencies,  however, 
we  can  never  exceed  a total  average  power  1.  5.  This  conclusion  can  be 
demonstrated  by  considering  a modulatio.  :tion  with  two  frequencies: 


M "(t) 


l + aj  sin  (Wj  t + (j^)  + a^  sin  t + 


(19) 


where  a^  and  a^  are  positive,  and  <t>^  and  ^ are  random  phases  greater 

II 

than  zero,  but  less  than  2tt.  The  requirement  that  M (t)  must  always  be 

8 

non-negative  means  (a^  + a^)  - 1.  Using  eq.  (16)  gives 


TS 


1 

7 


1+4-  (a.  + a 


/>  < 1 + 2-  <»,  + ^Z)Z  * 1.5 


(20) 


By  induction,  with  N frequencies,  P-j-g  can  never  increase  beyond  1.5. 

Since  the  dc  power  is  always  1.  0,  the  ac  power  will  never  exceed  0.  5. 

In  a multidither  COAT  servo  we  are  interested  only  in  speckle  modu- 
lations at  or  near  the  dither  frequencies,  because  only  they  can  significantly 

effect  the  convergence  process.  Thus  we  define  a "speckle  coefficient",  C , 

s 

as 


c 


s 


w.+A 

I 

w.  -Aw 


(21) 
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This  coefficient  is  the  sum  of  2N  integrals,  where  N is  the  number  of  dither 

frequencies  in  the  COAT  servo.  The  summation  is  carried  out  over  both 

positive  and  negative  values  of  the  dither  frequencies,  The  quantity  Aw 

is  the  open  loop,  unity-gain  bandwidth  of  each  COAT  servo  channel.  These 

quantities  are  shown  graphically  in  Fig.  23  for  the  positive  half  of  the  power 

spectrum  of  an  arbitrary  modulation  function.  The  shaded  areas  correspond 

to  the  integrals  in  eq.  (21),  while  the  total  area  under  the  curve  (including 

the  negative  frequencies)  corresponds  to  the  integral  in'eq.  (17). 

Physically,  C ' corresponds  to  that  portion  of  the  ac  modulation 

power  in  Mg(t)  which  is  at  or  near  enough  to  each  dither  frequency  such  that 

it  can  create  false  error  signals  which  compete  with  the  correct  dither  signals 

for  control  of  the  convergence  process.  Clearly,  the  sum  of  the  integrals 

in  eq.  (21)  will  always  be  less  than  or  equal  to  the  integral  in  eq.  (17).  In 

physical  terms,  this  is  the  obvious  statement  that  the  portion  of  the  ac 

speckle  modulation  power  at  or  near  the  dither  frequencies  is  always  less  than 

or  equal  to  the  total  ac  power.  Since  the  total  ac  power  is  less  than  or  equal 

to  0.  5,  then  C must  also  be  less  than  or  equal  to  0.  5.  Thus 
s 

C < -s/075  = 0.  707  . (22) 

s 

One  must  be  careful  not  to  confuse  the  power  in  the  speckle  modulation 

function  M (t)  (i.  e.  , eqs.  (16)  and  (17))  with  the  speckle  signal  power  cal- 
s 

culated  in  the  statistical  analysis.  The  former  does  not,  by  itself,  represent 

a physical  energy  rate.  It  is  just  a dimensionless  amplitude  function  which 

modulates  the  usual  "no-speckle"  intensity  at  the  receiver.  The  latter  refers 

to  an  actual  electric  signal  within  the  COAT  servo  which  falls  within  the 

2 

control  bandwidth.  It  is  directly  proportional  to  Cg  and  the  square  of  the 
mean  intensity  incident  on  the  receiver  as  shown  in  (eq . (26)). 
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The  generalization  of  the  statistical  analysis,  as  described  above 
results  in  new  expressions  for  the  servo  dither  signal,  Sp,  and  the  servo 
speckle  signal,  Sg: 


(23) 


Sc  = M KI 
S ac  m 


(24) 


Both  Sp  and  Sg  represent  physical  voltages,  with  the  following  definitions: 

K * photomultiplier  constant 

= mean  "no-speckle"  intensity  at  receiver 
Ip  35  dither  modulations  at  receiver 
M^c  * average  value  of  Mg(t)  = 1.0 
M = ac  part  of  M (t)  at  or  near  dither  frequencies. 

3.C  S 

The  rest  of  the  analysis  continues  in  the  same  manner  as  before.  We  find 
the  average  power  in  the  dither  signal,  P^,  and  the  average  power  in  the 
speckle  induced  signal,  Pg,  by  calculating  the  expectation  value  of  the  squares 
the  voltages  Sp  and  Sg.  The  results  of  these  calculations  are  given  in 
eqs.  (25)  to  (28). 


(25) 


» 


(26) 
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where 


and  where 

4*  = dither  amplitude  in  radians 

N = number  of  dither  channels 

a = the  angle  characterizing  the  convergence  level  in  the 
statistical  model^ 

J.(4j)  = the  ith  Bessel  function. 
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Figure  24  shows  the  analytically  predicted  maximum  and 
minimum  mean  convergence  levels  as  a function  of  the  speckle  coefficient 
Cfl.  As  shown  elsewhere,^  the  convergence  level  oscillates  between  two  well 
defined  levels.  It  should  be  pointed  out  that  these  analytical  predictions  rep- 
resent the  minimum  performance  that  can  be  expected  from  an  18-channel 
multidither  COAT  array  that  has  dither  amplitudes  of  0.349  radians  (±20°)  in 
each  channel.  There  are  a number  of  techniques  which  are  mentioned  later  in 
this  report,  and  which  are  still  under  investigation,  that  will  considerably 
improve  this  predicted  behavior. 

B.  Computer  Simulations 

During  the  past  quarter  we  received  from  the  General  Research 
Corporation,  three  computer-generated  time  functions  representing  receiver 
amplitude  modulations  due  to  speckle  through  an  arrangement  coordinated  by 
the  contract  monitor  Robert  Ogrodnik.  The  three  functions  were  generated 
by  computing  the  intensity  of  radiation  incident  on  a point  receiver  as  a 
function  of  time.  The  incident  radiation  was  modeled  to  correspond  to  the 
reflected  radiation  from  a spherical  target  (1  meter  radius)  illuminated 
uniformally  over  a 10  cm  diameter  spot  by  a coherent  10.6  pm  source.  Tbs 
receiver  is  co-polarized  with  the  transmitter  and  2 km  away  from  the  target. 
When  the  target  rotates,  a speckle  pattern  sweeps  by  the  receiver  causing 
amplitude  fluctuations  as  a function  of  time.  Thus  the  normal  receiver 
signal  containing  the  dither  modulations  is  multiplied  by  another  modulation 
function  representing  these  amplitude  fluctuations  as  shown  in  the  schematic 
in  Fig.  25. 

In  the  computer  simulations,  the  normal  receiver  signal  was  multi- 
plied by  one  of  the  three  time  functions  generated  at  GRC.  The  three  different 
functions  correspond  to  three  rotation  rates:  0.01  rad/sec,  2 rad/sec,  and 
10  rad/sec.  Figures  26  to  28  show  time  histories  and  Fourier  amplitude 
transforms  of  these  functions.  The  time  histories  have  been  normalized  by 

dividing  the  original  functions  by  5 so  that  M (t)  < 1;  the  Fourier  transforms 

s 

correspond  to  the  original  functions  as  generated  by  GKC.  The  dither  band 
used  in  the  computer  simulations  is  also  indicated  in  the  figures. 


60 


CONVERGENCE  LEVEL.  Iy 


Figure  25.  COAT  servo  simulation  showing  addition  of  "speckle"  noise  generator 
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Figure  26.  Computer  generated  speckle  modulation 

function  supplied  by  GRC  (£2  = 0.01  rad./ 
sec , ). 
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Figure  27.  Computer  generated  speckle  modulation 
function  supplied  by  GRC  ( ft  = 2 rad./ 
sec. ). 
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We  have  some  reservations,  which  we  have  communicated  to  GRC, 
about  how  realistic  these  modulations  really  are.  The  first  observation  is 
that  the  modulation  functions  shown  in  Fig.  27  and  2C  must  have  been  produced 
by  a speckle  pattern  with  a contrast  ratio  approaching  unity  because  of  the 
very  large  (essentially  100%)  modulation  depths.  Theoretical  and  experi- 
mental studies  done  by  N.  George,  3 together  with  experimental  evidence 
gathered  on  this  contract  and  discussed  in  Section  III,  indicate  that  such  high 
contrast  ratio  speckle  patterns  are  not  typically  observed  at  visible  wave- 
lengths. Since  targets  are  e''en  more  specular  at  IR  wavelengths,  it  is 
unlikely  that  such  patterns  will  be  typical  at  10.6  pm  either.  We  also  note 
that  the  modulation  functions  in  Figs.  27  and  28  repeat  with  a period  on  order 
of  a 0.  5 msr  c.  Ihis  is  clearly  unphysical  since  the  modulations  should  oe  a 
random  .unction  of  time. 

Another  reservation  we  have  communicated  to  GRC  concerns  th  ; use 
of  a point  receiver  (i.  e.  no  aperture  averaging)  in  the  computation  of  the 
modulation  function.  Realistically,  a finite— size  receiver  will  average  out 
some  of  the  intensity  fluctuations  and  thus  reduce  the  effect  of  hi<>h  contrast 
ratios.  We  are  still  in  the  process  of  resolving  these  questions  at  the  time 
of  this  writing,  but  we  anticipate  no  difficulty  in  establishing  a mutual  consen- 
sus between  ourselves  and  the  GRC  personnel  before  the  end  of  this  contract. 
The  results  will  be  documented  in  t\e  final  report. 

A speckle  coefficient  Cg,  was  calculated  for  each  of  three  GRC  cases 
in  order  to  compare  them  with  the  predictions  made  by  the  statistical 
analysis.  The  results  are  shown  in  Table  4.  The  large  uncertainty  in 
cases  2 and  3 arises  from  numerical  problems  in  computing  the  integrals  in 
eq.  (21)  for  the  spectrums  in  Figs.  27  and  28;  it  is  very  difficult  to  digitally 
compute  the  area  under  a curve  that  approaches  a delta  function.  A more 
continuous  spectrum  w*.,uld  be  easier  to  handle  and  more  realistic  as  well. 
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Table  4.  Speckle  Coefficients  of  the  GRC  Supplied 
Speckle  Modulation  Functions 


C 

s 

— 

Case  1. 

Q 

= 0. 01  rad/sec 

0.  04  (±0.  02) 

Case  2. 

n 

= 2.0  rad/ sec 

0.  6(±0.  05) 

Case  3. 

n 

= 10. 0 rad/sec 

0.  5 (±0.  10) 

Recall  that  C = 0.  707  is  the  absolute  worst  possible  value  consistent 
8 

with  conservation  of  energy.  This  corresponds  to  100%  of  the  ac  power  in  the 
speckle  modulation  functions  passing  through  all  the  COAT  servo  filters 
completely  unattenuated.  Note  that  cases  2 and  3 in  Table  4 have  values  of  C8 
that  approach  this  limit.  We  therefore  conclude  that  these  two  cases  are 
very  near  to  the  worst  possible  kinds  of  modulation  functions  — realistic  or 
otherwise. 

In  spite  of  these  reservations,  the  data  supplied  by  GRC  still  proved 
useful.  For  example,  we  know  that  the  divider  AGC  used  in  the  computer 
simulations  does  not  represent  the  actual  experimental  hardware  as  accurately 
as  we  would  like.  In  particular,  we  have  noti  'ed  that  in  simulation  runs  with 
very  large  amplitude,  narrow  band,  speckle  modulation  functions  centered 
on  one  dither  frequency,  there  was  a definite  lack  of  correlation  with  experi- 
mental data.  The  simulations  always  predicted  poorer  performance.  This 
was  attributed  partly  to  the  divider  AGC,  and  partly  to  the  signal  clipper  in 
the  servo  loop.  Therefore,  in  order  to  have  any  confidence  in  the  computer 
simulations  when  using  the  GRC  data,  we  first  had  to  investigate  the 
problems  in  the  computer  simulation.  This  investigation  has  shown  that  the 
simulation  AGC  creates  a significant  oecond  harmonic  signal  when  the 
speckle  modulations  approach  about  65%  modulation  depth.  We  have  also  been 
able  to  demonstrate  that  stronger  signal  clipping  in  the  servo  loop  makes  the 
simulations  correlate  much  better  with  experiment. 
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By  formulating  a better  model  for  the  AGC,  and  by  increasing  the 
clipping  in  the  servo  loop,  we  have  been  able  to  demonstrate  very  good 
agreement  between  computer -predicted  COAT  convergence  levels  and 
experimentally  observed  levels.  This  new  configuration  is  still  not  optimum, 
however,  since  it  predicts  generally  longer  convergence  times  than  those 
observed  without  these  changes.  This  investigation  will  therefore  continue 
into  the  next  quarter. 

With  the  above  reservations  in  mir.d,  we  now  present  the  results  of  the 
computer  simulations  using  the  speckle  modulation  functions  generated  by 
GRC.  Figures  29  to  31  show  time  histories  of  the  radiant  intensity  on  a 
target  glint,  normalized  to  unity  at  maximum  intensity,  for  the  case  where 
the  target  rotation  rates  are  0.01  rad/sec,  2.0  rad/sec,  and  10  rad/sec. 
These  cases  represent  the  best  results  achieved  so  far.  We  expect  the 
COAT  system  performance  to  improve  when  the  questions  about  the  AGC  and 
signal  clipping  raised  above  have  finally  beer  resolved,  and  when  a finite 
receiver  is  used  in  the  GRD  model. 

It  is  interesting  to  compare  the  predictions  of  the  statistical  analysis 

5 

developed  under  a previous  contract  (F30602 -75 -C-0001 ) and  generalized  in 
the  previous  section  of  this  report,  to  the  computer  simulations  just  dis- 
cussed. Using  the  values  of  C in  Table  4,  we  can  compare  the  analytical 

s 

convergence  levels  shown  in  Fig.  24  with  those  observed  in  Figs.  29  to  31. 

We  note  that  the  statistical  analysis  predicts  much  lower  convergence  levels 
fcr  cases  2 (2  rad/sec)  and  3 {10  rad/sec),  than  th^se  observed  in 
Figs.  30  and  31,  white  in  case  1 (0.01  rad/sec)  the  prediction  is  about  the 
same  as  the  observed  level  in  Fig.  29.  The  reason  for  the  agreement  or 
nonagreement  varies  depending  on  which  case  we  consider.  Case  3 for 
example,  has  two  narrow -band,  high -amplitude  frequency  modulations  within 
the  COAT  dither  band  (see  Fig.  28).  Since  the  statistical  ar  alysis  was  based 
on  a random  distribution  of  dither  modulations  at  or  near  all  the  dither  fre- 
quencies, one  would  not  expect  accurate  predictions  for  such  a spectrum. 
Cases  1 and  2 have  a more  continuous  spectrum.  Thus  we  would  expect 
better  correlation  between  prediction  and  observation.  Case  1 does  show  good 
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correlation,  while  case  2 does  not.  The  correlation  would  be  much  better  in 
case  2 if  signal  clipping  had  not  been  included  in  the  computer  serve  loop. 

To  demonstrate  this  conclusion,  Fig.  32  shows  a time  history  of 
case  2 (i.  e.  , using  the  same  speckle  modulation  function  as  Fig.  36),  except 
that  the  signal  clipping  is  left  out.  This  case  fails  to  converge  at  all,  which 
is  consistent  with  the  analytical  prediction.  Since  the  statistical  analysis 
does  not  include  the  effect  of  signal  clipping  on  the  power  ratio,  we  can 
always  expect  to  improve  upon  the  analytical  predictions  by  including  it  in  the 
servo. 
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Figure  32.  Simulation  with  GRC  modulation  function 
{SI  = 2.0  rad. /sec.)  (No  signal  clipping). 
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V. 


PJLANS  FOR  THE  REMAINDER  OF  THE  PROGRAM 


The  final  contract  quarter  will  be  devoted  to  experimental  measure- 
ments of  COAT  system  performance  with  the  speckle -characterized  targets. 
The  propagation  conditions  will  include  turbulence  and  thermal  blooming. 

Beam  convergence  level,  convergence  stability,  system  convergence  time, 
distortion  compensation  performance,  and  tracking  performance  will  be 
observed.  The  speckle  noise  spectrum  will  be  recorded  and  compared  with 
that  obtained  from  backscatter  measurements.  An  acousto-optic  speckle 
simulation  will  be  performed  for  the  purpose  of  assessing  the  minimum  modu- 
lation index  required  to  overload  a single  servo  channel. 

For  the  computer  simulation  studies,  we  will  continue  to  investigate 
techniques  which  improve  the  simulation  model  response  in  the  presence  of 
speckle  noise  and  also  improve  its  correspondence  with  the  real  system. 

New  speckle  data  to  be  generated  by  GRC  will  be  used  to  test  the  system 
response  and  convergence  parameters. 
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